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We report the first observation of an excited singly-charmed baryon SI* (ess) in the radiative 
decay fl°7, where the S2® baryon is reconstructed in the decays to the final states 0~ir + , S2~ty + tt°, 
S2~tt + tt~tt + , and S~ K~ tt + tv + . This analysis is performed using a dataset of 230.7 fb _1 collected 
by the BABAR detector at the PEP-II asymmetric-energy B Factory at the Stanford Linear Ac- 
celerator Center. The mass difference between the Q* and the S2® baryons is measured to be 
70.8 ± 1.0 (stat) ± 1.1 (syst) MeV/c 2 . We also measure the ratio of inclusive production cross 
sections of i7* and in e + e~ annihilation. 

PACS numbers: 13.30.Ce, 14.20.Lq 



The production of charm baryons is largely unexplored 
and provides an interesting environment to study the dy- 
namics of quark-gluon interactions. All singly-charmed 
baryons having zero orbital angular momentum have 
been discovered except for the J p = | + ess state, de- 
noted as SI*. A non-relativistic QCD effective field theory 
calculation predicts the difference between the mass of 
SI* (Ma* ) and the mass of Q° c (M n o), AM, to be between 
50 and 73 MeV/c 2 0. A lattice QCD calculation gives 
AM = 94 ± 10 MeV/c 2 0]. New quadratic baryon mass 
relations predict a mass of Mq* = 2767 ± 7MeV/c 2 
and several other predictions for Mq* exist around 2770 
MeV/c 2 [5-11], implying AM = 70 -"75 MeV/c 2 . 

Here we report the observation of an excited baryon SI* 
produced inclusively in e + e" — > SI* X processes, where X 
denotes the rest of the event. We measure the mass differ- 
ence, AM, and the ratio of the production cross section 
of e+e~ -> Q*X relative to e+e~ -> Q a c X. Throughout 
this paper, for any given mode, the corresponding charge 
conjugate reaction is also implied. 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric-energy e + e~ 
storage rings. The dataset corresponds to an integrated 
luminosity of 209.1 fb _1 collected at a center-of-mass 
(CM) energy of y/s = 10.58 GeV, near the peak of the 
T(4S) resonance, and 21.6 fb _1 collected approximately 
40 MeV below the Y(4S) mass. 

The BABAR detector is described elsewhere [l2| . 
Charged tracks are reconstructed with a five-layer, 
double-sided silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) with a helium-based gas mixture, 
placed in a 1.5-T uniform magnetic field produced by a 
superconducting solenoidal magnet. Kaons, pions and 
protons are identified using likelihood ratios calculated 
from the ionization energy loss (dE/dx) measurements 
in the SVT and DCH, and from the observed pattern of 
Cherenkov light in an internally reflecting ring imaging 
detector. Photons are identified as isolated electromag- 
netic showers in a CsI(Tl) electromagnetic calorimeter 
(EMC). Large samples of Monte Carlo (MC) simulated 
data are used for determination of signal detection effi- 
ciencies and for the optimization of the selection criteria. 
These are generated using JETSET [LJand the detector 
response is simulated with GEANT4[14T| . 

The SI* candidate is identified through its radiative 
decay, SI* — > L?°7, where the f2° is reconstructed ex- 



clusively in the following four decay modes, which are 
expected to provide the best signal-to-background ratio: 

• Q° c -> AK- (01) 

• f2° -> Sl-n+TT , Q~ -> AK- (02) 

• J?° -► L?-7r+7r-7r+, SI' -> AK~ (03) 

• J?° -> »-if-7r + 7r + , S~ -> Aw' (CI) 

The labels in parentheses to the right of each decay 
mode designate the four final states of the S7® decay. 

A A — > p7r~ candidate is reconstructed by identifying 
a proton track, combining it with an oppositely-charged 
track identified as a 7r~, and fitting the tracks to a com- 
mon vertex. Here and throughout this analysis, all re- 
constructed baryon candidates are required to have an 
acceptable x 2 from the vertex fit. The flight distance of 
each A candidate between its decay vertex and that of 
its parent (Sl~ or E~) is required to be greater than 0.30 
cm. The A — » pn~ signal is fitted using a sum of two 
Gaussian functions with a common mean. The signal re- 
gion is defined by |M p7r - -M A \ < 3.8 MeV/c 2 (» 2ct RM s), 
where Ma is the fitted peak position of the A and ctrms 
is defined by c^ms = /i <J i + H°2i where /i and fa are 
the fractions of the two Gaussian functions, and o\ and 
(72 are the two corresponding widths as obtained from 
the fit. The reconstructed A candidate is then combined 
with an identified K- (tt~) to form an Q- (E~ ) can- 
didate. The A and the K- (if-) tracks are fitted to a 
common vertex, and the flight distance of each L?~ or E~ 
candidate between its decay vertex and that of its parent 
(J7°) is required to be greater than 0.25 cm. Mass win- 
dows of \M AK - - M a - I < 5.2 MeV/c 2 (» 2ctrms) and 
\M A „- - M 3 - | < 6.0 MeV/c 2 (» 2a RMS ) are used to 
select Q- — > AK- and E~ — > An- candidates, respec- 
tively, where Mq- and Ms- represent the fitted peak 
positions of fl~ and Z~ . 

For the decay mode 02, the 7r° candidates are recon- 
structed by combining two photons. To enhance the ir° 
signal over combinatorial background, we require pho- 
tons to have a minimum energy of 80 MeV in the lab- 
oratory frame, to have a lateral shower shape consis- 
tent with that of a photon and to be well-separated 
from other tracks and clusters in the EMC. We require 
|M 77 - M^o] < 12.5 MeV/c 2 (2.5a), where M^o is the 
fitted peak position of the invariant mass of the two pho- 
tons. 
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For decays 01, 02, 03, the reconstructed Q~ is com- 
bined with a (ir + , ir + ir°, it + ir~it + ) to form an f2®, and 
fitted to a common vertex. For CI, the reconstructed S~ 
is combined with an identified K~ and two tt + tracks 
and fitted to a common vertex. The invariant mass 
of reconstructed candidates is required to lie within 
±2.5 ctrms of the central fitted value. The mass reso- 
lution is ctrms ~ 6 MeV/c 2 for 01, 03, and CI, and 
ctrms » 13 MeV/c 2 for 02. The resolution in 02 is dom- 
inated by the measurement of the photon energies from 
the 7t° decay. 

An Q* candidate is formed by combining a recon- 
structed f2® with a photon, applying the same photon 
selection requirements listed above for photons from ir° 
decay. For 02, it is required that the photon is not one 
of the 7T° daughters. 

Though eliminating most Q* baryons from B decays, 
the requirement that the scaled momentum of f2* candi- 
dates, (x p (/?*)), be greater than 0.5 significantly reduces 
combinatorial background from e + e~ —> qq (where q = u, 
d, s). The scaled momentum is defined as x p — P* /Pmax? 
where p* is the reconstructed momentum in the CM 
frame and p* max = -J 's/4 — M 2 , with M being the mass 
of the particle. 

Fig. [H shows the reconstructed invariant mass distri- 
butions of candidates with x p {t2®) > 0.5. Clear 
peaks indicating production of J?J? are visible in each of 
the modes represented in Fig. [TJ The invariant mass 
resolution is improved by 25% by using the variable 
M n - n + - M a - + Mf 2 DG , instead of M u - n+ , where M n - 
is the reconstructed mass of the f2~ and MqP is the 
world average mass of the f2~ [l[ . An unbinned extended 
maximum likelihood (ML) fit is performed to extract the 
signal yield. For each mode, a double Gaussian func- 
tion with a common mean is used to fit the signal and 
a first-order polynomial is used to model the combina- 
torial background. The mass resolution in each decay 
mode is obtained from a large sample of MC signal events 
reconstructed and processed in the same way as data. 
For the fits shown in Fig. [TJ the widths of the signal 
lineshapes are fixed to the values from MC simulation. 
The fit shown in Fig. [TJa) results in a raw (i.e. uncor- 
rected) yield of 156± 15 (stat) events and a mean mass of 
2693.3 ±0.6 (stat) MeV/c 2 . For the other three f2° decay 
modes the mean masses are fixed at 2693.3 MeV/c 2 , and 
a second-order polynomial is used to model the combina- 
torial background. The fitted raw yields are 92+ 2 g (stat), 
23±g° (stat) and 34±|| (stat) events for 02, 03 and CI 
decay modes, respectively. 

For Q* candidate selection, we require x p {Q*) > 0.5 
but make no direct cut on x p (f2®). The invariant mass 
distributions of Q* — > i?J?7 candidates are shown in 
Fig. [21 The invariant mass resolution is improved by 
w 40% by using the variable Mqo 7 - M n <> + Mf 2 P G , in- 
stead of Mqo 7 , where Mqo is the reconstructed mass 
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M a° - M \ + M ™° (GeV/c 2 ) 

FIG. 1: The invariant mass distributions of candidates 
reconstructed in the decay modes into (a) Q~ -n + , (b) 
n-7v + n°, (c) n-Tv + n-Tv + , and (d) S~K~w + w + . For all 
of these, we require x p (Q®) > 0.5. Here M n o is the re- 
constructed mass of the candidates, and denotes the 
daughter hyperon. The points with error bars represent the 
data, the dashed line represents the combinatorial background 
and the solid line the sum of signal and background. 



of the f2® and Mj^P G is the world average mass of the 
Q° (2697.5 MeV/c 2 ) [1]. A clear peak from Q* c -> J?° 7 
(J7j? — * i?~7r + ) production can be seen in Fig. EJa). The 
scaled sidebands which are also shown in Fig. [21 show 
no peak in the mass distribution. The distribution is fit- 
ted with the Crystal Ball function [l5[ to model the sig- 
nal and the product of a fourth-order polynomial and a 
two-body phase space function [l[ to model the combina- 
torial background. The signal shape parameters are fixed 
to the values found from MC simulation except for the 
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FIG. 2: The invariant mass distributions of Q* c — > i7°7 can- 
didates, with reconstructed in the decay modes (a) i? _ 7r + , 
(b) i7"7r + 7r°, (c) Q-^t^-^, (d) ~-if-7r + 7r + , and (e) for 
the combined decay modes (01, 02, 03 and CI). For all 
of these, we require x p (fl*) > 0.5. Here M n o 1 is the recon- 
structed mass of the ft* candidates, and M n o is the recon- 
structed mass of the • The points with error bars represent 
the data, the dashed line represents the combinatorial back- 
ground and the solid line the sum of signal and background. 
The shaded histograms represent the mass distribution ex- 
pected from the mass sideband of Q°- 



mean of the distribution. The invariant mass resolution 
is 4.0 MeV/c 2 . The fit results in AM = 69.9 ± 1.4 (stat) 
MeV/c 2 and a raw yield of 39tg° (stat) events. The fit 
is superimposed on Fig. [DJa). The signal observed for 
Q* — ► f2®j — ► j?~7r + ) corresponds to a significance 
of 4.2 standard deviations (a) including the systematic 
uncertainty on the observed yield. The significance is de- 
rived from y/ '2\n(L max / Lq) , where L max and Lq are the 
likelihoods for fits with and without a resonance peak 
component, respectively. The systematic uncertainty is 
discussed later. We use a similar fit procedure for 02, 03, 
and CI decay modes to extract the signal yields. For 03, 
Mqo is fixed to the value obtained from the process 01. 
The fits result in raw yields of 55j^g (stat), — 5±5 (stat), 
and 20±9 (stat) events for 02, 03, and CI, respectively. 

For all decay modes we determine the ratio of inclusive 
production cross sections, 

cr(e+e- -> Q*X,x p {Q* c ) > 0.5) 
~ a(e+e- -> f2°X, x p {{2° c ) > 0.5) ' 

where the scaled momentum of the Q* (f2®) is required 
to be greater than 0.5 in the numerator (denominator) 
cross section. We assume that B(f2* — > J?° 7) = 100%, 
and include 4?J? baryons coming from fl* decay as part of 
the denominator cross section, provided they satisfy the 
x p (S7®) requirement. The relative detection efficiencies 
(eo*/eQo) of the D* compared to fP c within these mo- 
mentum ranges are estimated from MC simulation and 
are listed in Table HI along with the results for the cross 
section ratios R. 

We combine 01, 02 , 03, and CI and perform a sin- 
gle ML fit. The fit results in AM = 70.8 ± 1.0 (stat) 
MeV/c 2 , a raw signal yield of 105 ±21 (stat) events, with 
a significance of 5.2er (including systematic uncertainty), 
and a ratio R = 1.01 ± 0.23 (stat). This proce- 
dure weights the individual decay modes by the observed 
number of f2® baryons in the data, and results in the 
minimum overall error on the combined value of R. The 
results are summarized in Table [I] 

Several sources of systematic uncertainty in the fitted 
signal yields are considered. The largest uncertainties 
arise from the fits to the mass spectra. These are esti- 
mated by repeating the fits, varying the fixed parameters 
of the fitted signal functions by ±1 standard deviation 
and varying the functional parametrization of the back- 
ground. The systematic uncertainty on the yield from 
the combined J?* modes is 6%. The systematic uncer- 
tainty on AM is dominated by the photon energy scale 
and is 1.5%. This is estimated from the distribution of 
reconstructed masses of low-energy neutral pions. The 
uncertainty in the fitting procedure leads to a systematic 
uncertainty of 11% on the ratio R, measured from the 
combined modes. There are also systematic uncertain- 
ties of 1.8% from the photon reconstruction efficiency, 
and 1.4% due to the limited MC sample size. The uncer- 
tainties from tracking, particle identification, selection 
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TABLE I: The mass difference, AM = Mn* — Mqo (MeV/c 2 ), the fitted signal yield, Y (events), the f2* signal significance, S 
(in a), the relative detection efficiency, en* /ej?o, and the ratio of inclusive production cross sections, R, as defined in the text. 
The first uncertainty is statistical, and the second is systematic. 



Decay mode 


AM (MeV/c 2 ) 


Y (Events) 


S(a) 




R 


OI 


69.9 ± 1.4 ± 1.0 


39I9 ± 6 


4.2 


0.35 


0.7llo;Jg ±0.11 


02 


71.8 ± 1.3 ± 1.1 


55t" ± 6 


3.4 


0.34 


1.76t°;S±0.21 


03 


69.9 (fixed) 


-5±5±1 




0.33 


-omtoli ±0.13 


CI 


69A -2.0 ± 10 


20 ± 9 ± 3 


2.0 


0.35 


1.70t^o ±0.34 


Combined 


70.8 ± 1.0 ± 1.1 


105 ± 21 ± 6 


5.2 


0.34 


1.01 ±0.23 ±0.11 



of intermediate hyperon candidates, daughter branch- 
ing fractions [l| and luminosity approximately cancel in 
the ratio, since the Q* analysis uses the same selection 
and data sample as the analysis. The sensitivity to 
fragmentation modeling is negligible. A possible addi- 
tional uncertainty arises from multiple candidates found 
in 10% of the events in the data, usually due to a com- 
mon hyperon combined with alternative particles from 
the rest of the event to form i?* candidates. These are 
uniformly distributed in Mq* and are hence absorbed 
into the background parametrization, with no evidence 
for multiple candidates peaking in mass. 

In summary, we report the first observation of an ex- 
cited singly-charmed baryon fi* (ess) decaying to J?° 
and a photon, with a significance of 5.2cr, and mea- 
sure the mass difference between Q* and to be 
AM = 70.8 ± 1.0 (stat)±l.l (syst) MeV/c 2 . This is 
consistent with the theoretical prediction in [2, 4-11] and 
below that described in Q ■ We also measure the ratio of 
inclusive production cross sections, 

R = 1.01 ± 0.23 (stat) ± 0.11 (syst). 
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